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THREE BALLISTIC CAMERA DATA REDUCTION METHODS 
APPLICABLE TO REENTRY EXPERIMENTS 

By John E. Hogge 
Langley Research Center 

SUMMARY 

Ballistic cameras were used by the Langley Research Center to photograph a reentry 
test object against a star background in order to determine accurately the position, 
velocity, and deceleration of the object as it penetrates the earth's atmosphere. Three 
data reduction methods which are applicable to the determination of the reentry trajec- 
tories from the photographic records are discussed. One reduction method applies when 
the reentry image on each photographic plate can be approximated with a straight line. 
Another method applies when the reentry image is curved on each photographic plate. A 
third technique applies when the shutter of each of several cameras at different stations 
is synchronized in time. 


INTRODUCTION 

Reentry experiments employing multistage rockets were performed from the NASA 
Wallops Station launch site at Wallops Island, Va., to investigate the physical phenomena 
which occur during reentry of high-speed objects into the earth's atmosphere. The first 
several stages lifted a velocity package (the reentry stages) up to a predetermined alti- 
tude; then the remaining stages were sequentially fired downward so that the last stage 
after burnout served as the reentry test object. Ionization and dissociation of body mate- 
rial and atmospheric gases produced a visible reentry trace against a star background. 
Ballistic camera exposures of the reentry trace against a star background were used to 
determine the position, velocity, and deceleration of the reentry object as it penetrates 
the atmosphere. 

A number of ballistic cameras were employed at the three NASA camera stations 
at Coquina Beach, N.C., Wallops Island, and Sandbridge, Va. At least one ballistic 
camera at each site was equipped with a time -controlled shutter to provide a programed 
sequence of exposures of the reentry trace for velocity determination. A camera net- 
work, consisting of one ballistic camera at each site, were synchronized in time — that is, 
each of the shutters opened and closed at the same instant of time. If all three cameras 
were successful in acquiring the reentry trace, then photogrammetric triangulation for the 



position of the reentry object as a function of time could be performed. When the reentry 
photographs were available only from nonsynchronized cameras, then other reduction 
techniques could be used. 

The purpose of this report is to summarize the methods and procedures that are 
available at the Langley Research Center for reducing raw photographic reentry data to 
final trajectory parameters. This reduction capability was established primarily for the 
Trailblazer I and II series of reentry experiments which are described in references 1 
and 2. 

Each of the data reduction methods discussed herein is written in the FORTRAN 
(FORmula TRANslation) n language for the IBM 7094 electronic data processing system. 


A 1 ,!* 1 ,* 1 

X 2 ,ju 2 ,^ 2 


a s> 5 s 


x s>y s 


A,B 

-Ap,Bp 

b,e 

K Bb’ K Be 


B b 

B e 

BB^Be 

^Bb’^Bb’^Bb 


SYMBOLS 

direction cosines referred to equinox 1 
direction cosines referred to equinox 2 
right ascension and declination of a star 

star image plate coordinates measured relative to an orthogonal 
system (x,y) on plate 

camera or station 

photographic plates from cameras A and B, respectively 
reentry image points on Bp 

unit vectors defining rays through projection center of camera B 
and end points b and e on Bp 

reentry point corresponding to image point b on Bp 

reentry point corresponding to image point e on Bp 

plane defined by camera B and the reentry 

direction cosines of 
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X Be ,M Be^Be 

L B 

^BP’^BP’^BP 

S A 

R BP 

X AB’ Y AB ,Z AB 

r ab 

m 

M 

X Am 

R Am 

w m 

W 

X Am’^Am’ ^Am 
X Am’ Y Am’ Z Am 

dm 


direction cosines of K Bg 

unit vector along X-, Y-, and Z-axis, respectively 

arc length of the reentry as seen from camera B, radians 

direction cosines of a normal to plane BB^B e defined by 
camera B and the reentry 

perpendicular distance from station A to plane BB^B e 

unit vector along the normal to plane BB^B e 

relative rectangular coordinates of station B with respect to 
station A 

distance between stations A and B 

reentry image point on Ap corresponding to beginning of a shutter 
break in data reduction method I 

reentry point corresponding to image point m 

unit vector defining a ray through station A and points m and M 

slant range from station A to reentry point M 

angle between and R^ m (see fig. 4) 

intersection of with plane BBpB e 

direction cosines of K^ m 

rectangular coordinates of reentry point M with respect to 
X,Y,Z centered at station A 

spatial length corresponding to mth shutter break on reentry image 
on Ap 
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X,Y,Z 

h' 

4>i 

a 

Ri 

SRi 


UT 

A i 

*o > % 
X ri’^r] ,p r] 

Wo 

Wc 

a c ,5 c 

° b >^0 

x o>y 0 

%s ,ri s 


astronomical equatorial coordinate system centered at station A 
geocentric latitude of ith camera station 
geodetic latitude of ith camera station 

semimajor axis of reference ellipsoid for earth, 6378.388 km 

radial distance to ith camera station from center of earth 

altitude of ith camera station above mean sea level 

sidereal time at ith station corresponding to instant of reentry 

sidereal time at Greenwich for zero hours universal time on night 
of reentry 

universal time of reentry 

west longitude of ith camera station 

standard coordinates of an object 

direction cosines of £-axis 

direction cosines of 77 -axis 

direction cosines of an object ray 

direction cosines of central ray 

right ascension and declination of central ray 

right ascension and declination of an object 

measured plate coordinates of an object 

standard coordinates of star image 
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Ax, Ay 
x,y 

x i^i 

U 

Q 

i^i 

X U> Y U ,Z U 

X q> Y q> Z q 

D 

X Aj’^Aj’ Z Aj 

X',Y',Z' 

X A’Y A> Z A 

c 

v y P 


difference between measured and computed values of star image 
plate coordinates 

orthogonal coordinate system in plane of photograph; all points of 
interest on plate are measured with respect to this system 

plate measurements on reentry image 

line connecting center of projection of camera A with some 
arbitrary point on the reentry 

line connecting center of projection of camera B with a particular 
point on the reentry 

direction cosines of U 

direction cosines of Q 

coordinates of an arbitrary point u on line U 

coordinates of an arbitrary point q on line Q 

distance between point u on line U and point q on line Q 

Cartesian coordinates of reentry point relative to station A and 
corresponding to jth image measurement on chopped plate B 

camera axis system in data reduction method III 

location of camera A with respect to fixed reference system 
X,Y, Z in data reduction method in 

focal length of camera 

coordinates of plate center in data reduction method in, measured 
relative to x,y system 
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direction cosines relating camera-fixed axes and space-fixed axes 
in data reduction method III 

initial approximation to reentry point 

error equations in least-squares process of data reduction 
method III 

angles relating camera-fixed axes to space -fixed axes in data 
reduction method III 

spatial distance along reentry trajectory 
time 
velocity 
deceleration 

Superscripts: 

h hours of time 

s seconds of time 

IDENTIFICATION OF STELLAR AND REENTRY EXPOSURES 

A sketch of a reentry, photographed against a star background, is shown in figure 1. 
This sketch represents a photograph which was taken with a stationary camera equipped 
with a time -controlled shutter - that is, a shutter programed to open and close at spec- 
ified time intervals. Hence, both the reentry and stellar images consist of trails which 
are broken into time -dependent segments. Programed exposures of the star field are 
made before and after reentry, with the post-reentry exposures made in the reverse 
time sequence. As can be observed from the figure, breaks in the star trails during the 
reentry period are not resolved; this is due to the rapid shutter rate which was programed 
for the reentry mode. The time associated with each break on the star trails is deter- 
mined from the initial time of camera opening and the subsequent shutter rate. The time 
corresponding to the occurrence of reentry is recorded at the camera site. Figure 2 is a 
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sketch of a typical photographic plate from a streaked camera — that is, from a camera 
that is not equipped with a time -controlled shutter. Note that both the reentry and star 
images are continuously trailed on the plate. 

From each photographic plate utilized in the reduction to final trajectory param- 
eters, a number of selected reference stars are chosen in the vicinity of the reentry 
images for the purpose of determining constants of the camera system. These constants 
appear in the transformation from photographic plate coordinates to standard coordinates 
relative to a perfect projection plane in data reduction methods I and II and are subse- 
quently used to correct for errors and camera system distortions in the data reduction. 
The reference stars are utilized in data reduction method in in the determination of 
camera orientation and other physical parameters of the camera system. The right 
ascension and declination of each reference star relative to the equinox of the beginning 
of the year nearest to the time of exposure are needed in the reduction methods. These 
quantities are determined in the following manner. The overall star field on the photo- 
graphic plate is identified with the aid of appropriate star charts. By visually comparing 
the plate with the star chart, approximate values of the celestial coordinates are read 
from the gridded star chart. With these approximations, the stars are identified in the 
Boss star catalog and the right ascension and declination relative to the epoch of the 
catalog (1950.0) are established. After these coordinates are corrected for proper 
motion, they are written in terms of direction cosines and transformed from the catalog 
epoch to the beginning of the year nearest to the time of exposure. According to refer- 
ence 3, this transformation becomes 


X 2 


X X Xy X z 


1 

M 2 

1! 

Y x Y y Y z 


M 1 

' "1 

I 


Z X Zy Z Z 


*4 


where and X^,/i^, are the direction cosines of the star relative to the 

epoch of the catalog and the year nearest to exposure time, respectively. In terms of 
the right ascension Q! s and declination 6 S of the star, 

X* = cos 6 S cos Q! s 
/!■*■ = cos 6 S sin a s 
= sin 6 S 
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The elements of the transformation matrix are defined and tabulated for the beginning of 
each year from 1900-1980 in reference 4, page 109. 

In addition to the celestial coordinates (right ascension and declination) of each 
reference star, the relative rectangular coordinates x,y of reference star images and 
points along the reentry traces, measured in the plane of the photographic plate, are 
required. On chopped plates, the beginning and end of each shutter break on the reentry 
trace are measured. Some convenient break in the star trails is measured with time 
corresponding to that particular break noted. On streaked plates, the beginning or end 
of the star exposure is measured. The procedure for measuring streaked reentry images 
varies with the technique of reduction and is therefore discussed with each reduction 
method. 

Precision comparators are employed at the Langley Research Center to perform 
photographic plate measurements." These machines have the capability of measuring 
relative distance of points on the plate to within ±5 microns over the entire travel of the 
engine ways. 

DATA REDUCTION METHOD USING STRAIGHT-LINE REENTRY IMAGES 

The method which uses a straight-line approximation to the reentry image on each 
plate (herein referred to as data reduction method I) is essentially the method discussed 
in reference 3. This analysis technique was developed for the reduction of meteor trails 
observed simultaneously from two camera sites. In the solution for the spatial position 
of each shutter break, the assumption is made that the reentry is part of a great circle 
when projected onto the celestial sphere. If this assumption is valid, the reduction of 
photographic plate coordinates to final position, velocity, and deceleration can be per- 
formed with nonsynchronized pairs of photographs. Furthermore, a complete reduction 
can be made when only one of the photographic images is chopped. These features have 
led to extensive use of this analysis technique for the reduction of trajectory parameters 
from photographic observations. A general description of the method as it was applied 
to the reentry problem is given herein. Much of the details of the technique is omitted 
since a complete discussion can be found in reference 3. However, a derivation of some 
of the working equations is given in the appendixes herein to present the geometrical con- 
nection of the camera network parameters to these equations. 

To illustrate the technique, it is assumed that one of the photographic plates is not 
chopped and is arbitrarily chosen to be plate B in figure 3. Furthermore, it is assumed 
that the reentry image forms a straight segment on plate B and that reentry image breaks 
on plate A can be connected with a straight line. In figure 4, A and B denote the 
cameras from which the plates A p and B p were obtained. 
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To establish the plane BBjjB e through the center of projection of camera B 
and the reentry, two points on the reentry image are selected. One point should be near 
the beginning of the reentry image and the other point should be near the end as indicated 
by b and e, respectively, in figure 3. Unit vectors Kg^ and Kg e (fig. 4) define 
rays through the projection center and the end points b and e on the photographic 
plate. These rays pass through the reentry points Bfc and B e . The unit vectors can 
be written 


%b " iX Bb + ^Bb + ^Bb 

and 


K Be " 1 X Be + 


j^Be + kv Be 


where i, j, and k are unit vectors along the X, Y, and Z axes and Xg^, p g^, j/g^ 
and ^B e >Mg e >^g e are the direction cosines of Kg^ and Kg e , respectively. These 
direction cosines are relative to the reference coordinate system which at this point is 
defined only to be centered at camera A. 

Direction cosines of a normal to the plane BB^B e , defined by station B and the 
reentry points B^ and B e , are determined from Kg^ and Kg e to be 


^Bb^Be ~ ^Bb^Be 
BP sin Lg 


_ ^Bb X Be - X Bb^Be 
M BP - sin Lb 


_ X Bb M Be “ ^Bb X Be 

^BP “ sin L g 


( 1 ) 


where Lg is the arc length of the reentry in radians, as seen from station B, and is 
determined from observing that 


V BP 


M B p 2 + 


v 2—1 
^BP 1 
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The perpendicular distance from camera A to the plane BB^Bg is given by 



R AB _ X AB X BP + Y AB M BP + Z AB^BP 


where 


P BP " iX BP + ^ m BP + kt, BP 
R AB = lX AB + " y ab + kZ AB 

and ^AB ,Y AB’ Z AB are ^ ie re l a ti ye rectangular coordinates of camera B with respect 
to camera A. 

The unit vector K Am (fig. 4) defines a ray through the center of projection of 
camera A and some measured shutter break m on Ap. In terms of direction 
cosines, 


X Am = Tx Am + ^Am + ^Am 


The intersection of this ray with the plane BB^B e is the point on the reentry M whose 
image on A p is m. The coordinates of this point X Arn ,Y ATr| ,Z Air| are determined 
as follows. 

The slant range R Am from station A to point M is determined from the right 
triangle AWM in figure 4 by the equation 


R Am = S A sec w m 


( 2 ) 


where W is the intersection of with the plane BB^Bg and o> m is the angle 
between and R^^. From figure 4, 


P BP * X Am " cos - X BP X Am + ^BP^Am + ^BP^Am 


and hence 


_ X AB X BP + Y AB m BP + Z AB^BP 
Am X BP X Am + ^BP^Am + ^BP^Am 


( 3 ) 
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The components of R Arn along the reference axis are 

X Am = R Am X Am 
Y Am = R Am' J Am > 
Z Am = R Am P Am 


( 4 ) 


The application of equation (4) to the beginning and end of each shutter break on the 
photographic plate results in the X,Y, Z coordinates of the corresponding points in 
space relative to camera A. The spatial length corresponding to a particular shutter 
break on the plate is 


“m » [(JWi - X^) 2 ♦ (Y Am+1 - Y Am ) 2 + ( ZAm+1 - Z^) 2 ] V2 (5) 

where subscripts m and m+1 denote the beginning and the end of a break, 
respectively. 

In order to apply equations (1) to (5), direction cosines of points along the reentry 
trajectory as seen from both camera sites and the relative coordinates of station B 
with respect to station A must be determined. The direction cosines are determined 
from measurements on the photographic plates and the relative station coordinates are 
found from a knowledge of the station location (longitude, latitude, and altitude above sea 
level) and the time at which the reentry occurred. A discussion of each of these compu- 
tations follows. 

The axis system shown in figure 4 is an astronomical equatorial system centered 
at station A. The X-axis points toward the equinox of date and the Y-axis lies 90° east 
in the plane parallel to the equator. The Z-axis is chosen to form a right-hand coordi- 
nate system. Hence, the coordinates of station B with respect to station A are func- 
tions of time. To calculate these coordinates, certain auxiliary computations are neces- 
sary. These computations include the earth's local radius at each camera site, 
conversion of geodetic latitude to geocentric latitude, and the sidereal time at each 
camera station corresponding to the time of reentry. 

For the reduction of reentry photographic observations, the International Ellipsoid 
is used to represent the mean sea level of the earth. Based on this reference model, the 
geocentric latitude of and the radial distance from the center of the earth to the ith cam- 
era station are determined from the following equations from reference 5 (p. 489): 


11 


0/ = - 11* 3576635 sin 2^ + lVl731 sin 40 t - 0V0026 sin 6^ 


( 6 ) 


and 

R A = a^O. 998320047 + 0.1683494 x 10 -2 cos 2^ - 0.3549 X 10“ 5 cos 4^ 

+ 0.8 x 10“ 8 cos Q<Pij + 6Ri (7) 

where 

<p^ geodetic latitude of ith camera 

geocentric latitude of ith camera 

a semimajor axis of reference ellipsoid, 6378.388 km 

R^ radial distance from earth's center to ith camera station 

6Ri altitude of ith camera station above mean sea level 

The sidereal time at the ith station 0 ^ corresponding to the instant of reentry is 
computed from 


01=00+ (UT) - A i + A0 S 


where 

0q sidereal time at Greenwich for 0* 1 universal time on night of reentry 

UT universal time of reentry 

Aj west longitude of ith station 

A 0 s = 236.555(uTh/24 h ) 

s,h seconds and hours of time, respectively 

With these parameters determined, the relative rectangular coordinates of station B 
with respect to station A at the instant of reentry are computed from the following 
equations which are developed in appendix A: 
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X AB = |^B COS *B cos ( A B ' A A) “ R A COS ^a] cos 9 A" |^B COS *B sin ( A B " A A)] Sin 6 A 
Y AB = [*B cos *B sin ( A B - A A)] C0S °A + cos *B C0S ( A B " A A) " R A COS ^ sto 0 A 


> ( 8 ) 


Z AB ~ R B sin ^B ' R A sin *A 


The distance between the two camera stations is given by 

R AB - ( X AB 2 + Y AB 2 + Z AB 2 ) 172 

To calculate the direction cosines of points along the reentry trajectory, by utilizing 
the measured plate coordinates, the method of standard coordinates (see appendix B) is 
used. This method assumes that the photographic plate is a perfect projection of the 
object space; that is, points on the celestial sphere are projected through its center 
(projection center of the camera) onto a plane surface (the photographic plate). 

If the right ascension and declination of the central projection ray (i.e., the ray 
normal to the projection plane) and an object in space are known, the standard coordi- 
nates of the object (appendix B) are given by 


where 


£ 

° *o x c + M oMc + v o v c 

_ Vt? + Wy + 'Vi? 
° ^o^c + ^o^c + ^o^c 


■'v 


> 


( 9 ) 


X 0 = cos 6 0 cos Oq 
ju . 0 = cos 6 0 sin c/q 
v Q = sin 6 0 

X^ = -sin S c cos 
= -sin 6 C sin 

v tl = cos 6 c 


x c = COS 6 C COS C2q 

H c = cos 6 C sin q c 
v c = sin 6 C 

X^ = -sin 
/jL| = cos 6 C 
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and Qq, 6q and a^,6 0 denote the right ascension and declination of the central ray and 
the object, respectively. 

The £ and 77 coordinates are relative to the perfect projection plane. If the 
central ray is taken to be the plane center, then for each object whose right ascension 
and declination are known the corresponding standard coordinates can be computed from 
equations (9). Conversely, if the standard coordinates of an object are known, the direc- 
tion of the object in space can be determined by the inverse relationships 


^-c + ^77^0 + ^£^0 

\A + * 0 2 + 


= 


Mo “ 


v - 




Jl + £ q 2 + ?] 0 2 


V + V V 
C 7) 7 0 


i 1 * «o 2 + 


( 10 ) 


where X Q , ijl 0 , and v Q denote the direction cosines of the object with respect to an 
astronomical equatorial coordinate system centered at camera A. Hence, the direction 
of points along the reentry trajectory is determined once the corresponding standard 
coordinates are computed. 

It is assumed that the measured plate coordinates of an object x 0 ,y 0 are linearly 
related to the standard coordinates Z 0 ,V 0 by 


x o = ^o + Mo + c x 
= ay^o + ^^o + c y 


(ID 


where a x , bj,, c x , ay, by, and Cy are plate constants to be determined. Approxi- 
mately 10 reference stars, whose right ascension and declination are known, are selected 
in the vicinity of the reentry image on the photographic plate. The star coordinates are 
measured as described in reference 3. Therefore, both the measured coordinates x s ,y s 
and the standard coordinates % S ,V S are available since the right ascension and declina- 
tion for each reference star is known. With this information, the plate constants 
appearing in equations (11) are determined by the method of least squares. Once these 
plate constants are established, the inverse relationships are used to determine the 
standard coordinates and, hence, the direction cosines of any point of interest on the 
photographic plate from the comparator measurements. Therefore, the standard 
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coordinates of points along the reentry trajectory are determined from the photographic 
image measurements x 0 ,y 0 and the inverse of equations (11) 


where 


and 


«o = V o 
% = Vo 


+ Vo + % 
+ Vo + c 


V 


_ b y 

V> — 

-bx 

„ _ -Vx + b x c y 

c xy 

V 

c xy 

* c X y 

= 

t> = 

% 

ay C X - ^xCy 

V Cxy 

V 

c xy 

V C x y 


c xy = 

:a x b y ■ 

■ 3yb x 


( 12 ) 


Before the standard coordinates of reentry points are computed, the measurements 
on the reentry image are corrected as described in reference 3 for optical and emulsion 
distortions as well as for incorrect identification of the plate center, differential aberra- 
tion, differential refraction, and errors in the comparators. These corrections are per- 
formed as follows: 

With the plate constants determined, standard coordinates of reference 
stars are substituted into equations (11) and the corresponding plate coordinates 
are computed. The computed coordinates are compared with the measured 
coordinates and the differences Ax and Ay are used to correct the plate 
measurements of the reentry image. For the ballistic cameras employed in 
the optical coverage, the plate residuals Ax and Ay were on the order of a 
few microns. The plate residuals Ax and Ay are plotted separately 
against x along the reentry image and the points are connected with smooth 
curves. Corrections are read from these graphs and applied to the reentry 
image measurements. The standard coordinates are then determined from the 
corrected measurements through the application of equations (11), and then the 
direction cosines of the reentry points are computed from equations (10). 
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DATA REDUCTION METHOD USING CURVED REENTRY IMAGES 


Data reduction method I has the advantage that, if the reentry images on the plate 
can be approximated with a straight line, a complete reduction is possible when one 
chopped plate is used in conjunction with a streaked plate from a different camera site. 
Also, a complete reduction can be made when two nonsynchronized chopped plates are 
used. A modification to data reduction method I, applicable to the reduction of non- 
synchronized pairs of photographs where the reentry images show definite curvature, 
was made. A description of this technique (data reduction method n) with some of the 
working equations developed in appendix C follows. 

Sketches of a streaked photograph and a chopped photograph of a curved reentry 
are shown in figure 5. On the chopped plate, measurements of the beginning and end of 
each shutter break are made as data reduction method I. Reentry image measurements 
on the streaked plate are made at 250-micron intervals along the image. If the curva- 
ture of the image is smooth, a third-order polynomial is fitted through the measure- 
ments by the method of least squares; hence, for any xj along the image 

Yi = 0 O + h x i + ^2 x i 2 + ^3 x i 3 (* 3 ) 


where the constants /3 q, j3j, /3g which appear in equation (13) are determined from 
the least -squares fit. 

As described in data reduction method I, if the corrected plate coordinates xj 
and y^ are given, the corresponding direction cosines can be determined. By evalu- 
ating equation (13), at a few micron increments in the x-direction, a set of direction 
cosines corresponding to each increment can be tabulated. 


In figure 6, A denotes the camera corresponding to the streaked photograph. 

Let the line U emanating from camera A be defined by direction cosines iq 

obtained from some arbitrary point x^y^ on the curved streak on plate A. The line Q 
beginning at camera B and extending in space has direction cosines A j , ju. j , , which are 
determined from some measured shutter break Xj,yj on the chopped plate. 

The distance between lines U and Q is 


D =\/( X u-Xqf + 


u - ^( X q - X AB) - Y AB 




AB 


)- 


- 1 2 


J AB 


(14) 


As derived in appendix C, equation (14) yields the minimum value of D when 
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where 


X 


u 


_ x i( B i ~ C j A i.i) 


1 - A, 


13 


x q = 


_ X j( B i A ij ~ C j) 


+ X 


1 - A- 


AB 


i] 


B i " X i X AB + M i Y AB + p i Z AB 
C j = X j X AB + ^j Y AB + t ’j Z AB 
A ij = *i x j + MjMj + 

X^g, Y^g, are the relative station coordinates, and the subscripts i and j refer 
to plates A and B, respectively. 

Among all the sets of direction cosines computed by evaluating the poly- 

nomial fit to the curved streaked reentry, one set gives D a relative minimum value. 
This relative minimum value is determined by evaluating D for every set of direction 
cosines and choosing the smallest value. It is assumed that the set which 

gives rise to the relative minimum distance D defines a line which passes through the 
point on the reentry trajectory corresponding to the measured shutter break on the 
chopped plate. Hence, a triangulation is performed with the two rays and the 
vector R^g, which is shown in figure 6 and described in data reduction method I. This 
operation results in the spatial coordinates, relative to camera A, of the reentry point 
corresponding to the shutter break. The rectangular coordinates of the reentry point are 
given by 


v _ x i( B i “ c j A ij) 

A i ' 1 . 4. .2 




V 


Y Ar£r x Ai 


z Ar^r x Ai 


( 15 ) 
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and the slant range becomes 



C j A i.i 


A. 

A 1J 


These expressions for the spatial coordinates are evaluated by using the set of direction 
cosines A^/jq, that gave D a relative minimum value. By the application of this 
analysis technique to the beginning and the end of each shutter break on the chopped 
plate, the spatial components of each of these points on the reentry trajectory can be 
determined and, hence, the distance between points can be computed. 


DATA REDUCTION METHOD USING SYNCHRONOUS PHOTOGRAPHS 


For synchronized photographs, the data reduction method (herein referred to as 
data reduction method III) described in reference 6 is applicable. Since the details of 
the theory are developed therein, only a brief discussion of this technique with the working 
equations is given in this report. The notation for the working equations of reference 6 
presented herein is that of the present report. 

Figure 7 illustrates the connection between the image space (photograph) and the 
object space (reentry). The measured image coordinates x,y are related to the spatial 
coordinates X,Y,Z of the corresponding point on the reentry by 


x 1 (x-x a ) + m 1 (y-y a ) + -' 1 (z-z a ) 


x 3 (x - X A ) + 4 3 (y - Y A ) + * s (z - Z A )_ 



> 

x 2 (x-x a ).m 2 ( y -y a ) + ^( z -Z a ) 


h( x - *a) + " 3 ( y - Y a) + v i( z - z a) 

J 


(16) 


Two such equations, result from each camera in the network. The direction cosines 




’> v 3 


relate the camera-fixed axes X* ,Y',Z' to the fixed reference frame 


X,Y,Z which in this data reduction method is a local vertical system with origin at some 
desired latitude and longitude. The X-axis points north, the Y-axis lies 90° east in the 
tangent plane, and the Z^-axis lies along the local vertical. These direction cosines are 
functions of three angles Oj , c^, and cto which are the azimuth, elevation, and roll of 

the camera, respectively. These camera orientation parameters are referred to in ref- 
erence 6 as rotational elements of exterior orientation. The interior camera elements 
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Xp,yp and c are the coordinates of the photographic plate center and the focal length 
of the camera, respectively. The plate center coordinates Xp,y p are measured with 
respect to an axis system on the photographic plate. These six camera orientation con- 
stants are determined from a distribution of reference stars on the photographic plate. 
The technique for determining these constants by utilizing the reference stars is dis- 
cussed subsequently. From a knowledge of the latitude, longitude, and altitude above 
mean sea level of station A, X A ,Y A ,Z A with respect to the reference frame can be 
determined. 

Equations (16) provide two conditions on three unknowns - the X,Y,Z coordinates 
of the reentry point corresponding to the measured coordinates x,y. If the orientation 
constants and the position of each camera in the network are known, then three cameras 
in the network will provide six conditions on the coordinates of the reentry point. Hence, 
an overdetermination of the solution for these coordinates is possible. The method of 
differential correction (ref. 7, pp. 309-311), utilizing the least-squares scheme, is 
employed to solve these equations. To apply this technique, the equations are linearized 
by keeping only first-order terms in a Taylor's series expansion of equations (16) about 
an initial approximation to the coordinates of the reentry point. The error equations 
whose squares are minimized in the least-squares solution are 


e x = f( X °’ Y 0’ Z o) 
e y = g(Xo,Y 0 ,Zo) 


+-2L 

ay , af 
AX h 

ay + ^ 

ax 

8Y 

az 


o 

O 


1 AX + 

AY + 

ax 

o 

CD 

k! 

o 92 


AZ - x + Xp 
AZ - y + y p 




(17) 


where 


f 


g 


X 1 

( x - X a) 

+ 

( y - y a: 

+ u 1 

N 

it 

CO 

* 1 

^x - x ^ 
X X A, 

+ "s( Y - y a) 

1 + ^3 

( z - Z a). 


* 2 ( x - x a) 

+ p 2 

( y - 

+ v 2 

^isT 

"^1 

1 V 

CO 

1 

t> 

*»*3( Y - Y a) 

| + 

( z - z a)_ 


AX = X - X Q 
AY = Y - Y q 
AZ = Z - Z Q 
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X 0 ,Y 0 ,Z 0 is the initial approximation to the reentry point, and the o subscript to the 
partial derivatives imply that the partials are evaluated by using the initial approxima- 
tion Xo,Y 0 ,Z 0 . 

For one object point, two error equations e x and Cy exist for each camera. 

The most probable position of the corresponding reentry point is the one for which the 
sum of the square of the errors (e's) is minimized. 

Through the least-squares process AX, AY, and AZ are determined from the 
initial approximation. These quantities are added to the initial approximation to give a 
new approximation to X,Y,Z. Substitution of this new approximation into equations (16) 
shows the extent to which these equations are satisfied. If equations (16) are not satisfied 
within some prescribed limit, the computed values of X, Y, and Z are used as a new 
approximation and the process is repeated until convergence is obtained. The values for 
the variables which result in convergence are the most probable position coordinates of 
the reentry point. This technique is applied to the beginning and end of each shutter 
break on the photographic plate. Hence, the spatial coordinates of points along the reentry 
trajectory corresponding to each shutter break are determined. The spatial distance 
corresponding to each shutter break on the plate can now be computed. 

Equations (16) can be written in directional form by dividing both numerator and 
denominator of each equation by the magnitude of the slant range vector from the center 
of projection to the reentry point. The following equations result: 


x - x p = c| 


OijA + juj/x + 
K\ + p 3 p + 


y - y p = c i 


+ / u 2 M + 


(18) 


where A, p, and v are the direction cosines of the object relative to the fixed-axis 
system and are expressed, respectively, as 


X - X, 


A = 


R* 


Y - Yi 


M = 


v = 


R A 
Z - Zj 
R. 
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with 



(x - x A ) 2 + (y - y a ) 2 + (z - Z A ) 



In contrast to equations (16), equations (18) show no dependence on the absolute position 
of the object. 


For each reference star on the photographic plate, both the measured image coordi- 
nates x s ,y s and the corresponding direction cosines X s ,/i s ,^ s are known. These 
quantities are determined from the plate measurements and the right ascension and 
declination of each reference star as determined from the star catalogs and charts. From 
a knowledge of the station location and the time of star exposure,' the right ascension and 
declination of each reference star is transformed into the direction cosine relative to the 
fixed- axis system X,Y,Z. Hence, for n reference stars chosen on the photographic 
plate, 2n equations in six unknowns Xp, yp, c, Oj, Og, cr<j exist for each camera. Approx- 
imately 15 reference stars which are well distributed over each photographic plate are 
used to determine the elements of camera orientation for each camera. The method of 
solution is identical to that previously described for the position coordinates. The initial 
approximations to the azimuth, elevation, and roll of the camera a i> a 2 ,cr 3 are obtained 
from the camera dial readings, and the focal length c is approximated with that supplied 
by the camera manufacturer. The plate center coordinates Xp,yp are approximated by 
their values as computed from measurements of the fiducial marks on the photographic 
plate. The magnitudes of <7j, Og, °3> x p> Yp> an< * c which satisfy equations (18) to 
some prescribed limit are the most probable values for these camera constants. The func- 
tional relationships between o^jOgjOg and Xj,Xg, • • *>^3 are derived in reference 6. 


VELOCITY AND DECELERATION 


The spatial distance corresponding to the length of each shutter break on the plate is 
determined by the application of one of the foregoing methods. Elapsed time from the 
beginning to the end of a particular shutter break is determined when the period of revo- 
lution of the shutter and shutter configuration are known. After the distances along the 
reentry trajectory have been properly weighted (in a manner discussed in ref. 3) to 
account for image quality and readability on the plate, they are fitted by the least-squares 
method to the following equation which is found in reference 3: 


d = aj 


+ bjt 


+ Cje 


kjt 


(19) 


In equation (19), d and t denote the spatial distance along the reentry trajectory and 
the time, respectively. The constants aj,bj,Cj,kj are determined from the least- 
squares process. 
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Once the constants are determined, the velocity and deceleration equations result 
immediately by differentiation; thus 


kit 

V = bj + Cjkje 

• 9 kit 

V = c-jk^e 1 




( 20 ) 


J 

Hence, the velocity and deceleration at any time along the reentry trajectory are computed 
by evaluating equations (20) at the desired time. 


As described in reference 3, equation (19) is an empirical approximation to the 
true variation of distance with time along the reentry trajectory. For long reentry 
images, a better representation of distance as a function of time is obtained when the 
image is segmented and equation (19) is applied to each segment. 


CONCLUDING REMARKS 


The three data reduction methods discussed in this report were assembled for 
application to Langley reentry experiments in which photographic coverage was employed 
for position and velocity determination. The data reduction method using straight-line 
approximations to the reentry images (data reduction method I) has been applied exten- 
sively to the Trailblazer I and II series of reentry experiments discussed in NASA 
Technical Notes D-2189 and D-1866 because most of the reentry images obtained were 
approximately straight-line segments on the plates and in most of the experiments only 
nonsynchronous photographs were available. The data reduction method utilizing a 
third-order polynomial approximation to the reentry images (data reduction method II) 
has been applied to several pairs of plates from the Trailblazer experiments and the 
results indicate that this method of reduction should be more suitable than data reduction 
method I where the reentry image shows a definite curvature on the film. However, 
limited data requiring data reduction method II have been available and a full evaluation of 
this technique has not been possible. Although the data reduction method using synchro- 
nous photographs (data reduction method HI) has been checked out with hypothetical data 
and data from other projects, suitable data from the Trailblazer reentry experiments 
have not been available for the application of this technique. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., March 28, 1967, 
125-23-02-04-23. 
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APPENDIX A 


RELATIVE STATION COORDINATES AT INSTANT OF REENTRY 

In figure 8, R A and Rg denote the radii from the center of the earth of the two 
camera stations A and B, respectively. The respective geocentric latitudes are rep- 
resented by <p A and <P-q' . At the instant of reentry, let $ A be the sidereal time at 
station A. The X-axis points toward the mean equinox of date and the Y-axis lies 90° 
east in the plane of the equator. The Z-axis is chosen to form a right-hand orthogonal 
system. 

From figure 8, 


X A = R a cos <p A cos 0 A 
Y A = R A cos sin 9 A 


Z A = r A sin ^A 


X B = R B cos ^B* cos I 

: e A + y 

Yg = Rg cos <p-Q sin( 

e A + r) 

Z B = R B sin ^B 



where y is the difference in the west longitudes of each site or 

y = A b - A A 

Therefore, the coordinates of camera B with respect to camera A (eqs. (8)) are 


X AB = X B 

y ab = y b 

Z AB = Z B 
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or 

X AB = ( R B cos 0 B ' COS y - r a cos 0 A '^cos 0 A - ^Rg cos 0 B ' sin y^sin 6 A 
y ab = ( r b cos 0g’ sin y^cos 0 A + ^Rg cos $ g' cos y - r a cos 0 A '^sin 0 A 

Z AB = R B sin ~ R A sin ^A* 
and the distance between the camera sites is given by 

R AB - ( X AB 2 + Y AB 2 + Z AB 2 ) 1/2 
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APPENDIX B 


STANDARD PLATE COORDINATES 

The plane tangent to the celestial sphere at the point where the optical axis of the 
camera intersects the sphere is represented by the hatched area in figure 9. The photo- 
graphic plate is parallel to this tangent plane. (See ref. 8.) In figure 9, let C denote 
the center of the unit sphere, C' the point where the optical axis (central ray) intersects 
the sphere, O an object in the vicinity of C', O' the projection of the object on the 
tangent plane, and N the north celestial pole. The components of CO' along the 
£- and rj-axes are defined to be the standard coordinates of that particular object. The 
77 -axis increases in a northerly direction corresponding to increasing declination. The 
£j-axis is the projection of the parallel of declination through the tangent point and 
increases in an easterly direction corresponding to increasing right ascension. The 
origin of the standard coordinate system is the point of tangency, and the 77 -axis is the 
projection of a meridian through the tangent point onto the projection plane. 

From figure 9, 


tan^=^-=CO’ (21) 

where p is the angle between the optical axis and the object ray. Also, from the 
figure, 


£ 0 = C'O' sin Q = sin 0 tan p 
rj Q = C'O* cos 0 = cos 6 tan p 


( 22 ) 


If the direction cosines in an astronomical equatorial coordinate system of the object 
ray CO and the optical axis CC' are denoted by A 0 ,jli 0 ,v 0 and A c ,jn c ,^ c , respec- 
tively, then 


cos p = A 0 A C + p Q p c + u Q v c (23) 

where 


A 0 = COS S 0 COS 

A c = cos 6 C cos (% 

]U 0 = cos 6 0 Sin 0 ^ 

p c = cos 5 C sin 0 ^ 

v Q = sin 6 0 

u Q = sin 6 C 
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and Oq, 5 q and c^,6 c represent the right ascension and declination of the object ray 
and optical axis, respectively. 

By use of equation (23), equations (22) become 


t _ sin 6 sin \p 
4 ° ~ x o x c + M 0 M C + Wc 

cos 8 sin xU 
= 






0 x o v * * * * x c + ^c + ^c 
From the spherical triangle OC'N in figure 9, 

sin i p cos 8 = sin 6 0 cos 6 C - cos 6 0 sin 6 C cos^o^ - afc) 
sin i// sin 8- cos 6 0 sin(a 0 - 
From figure 10, the direction cosines of the £-axis are 


= cos 


jU^ = cos 


0 cos^, + = -sin oIq 

0 sin^, + |j 


cos 


v ^ = sin 0 = 0 


(24) 


(25) 


(26) 


and those of the rj -axis are 


X^ = cos [5 C + Ijcos Oq = -sin 6 C cos a Q 

jXq = cos ^6 C + l^sin = -sin 6 C sin \ 

v v = sin ( 5 c + |) = cos 6 c 

Substituting equations (26) and (27) along with expressions for \ 0 ,ii 0 ,v Q 
tions (25) yields 


(27) 


into equa- 
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sin ^ cos B - 
sin xj/ sin 0 = X^X 0 + 

Hence, substitution of equations (28) into equations (24) gives 

XqX^ + 

\>^c + ^o ^c + v o v c 

^ 0^77 + ^o^T) + Wv 
T] °' X Q X C + [£q (1q + v Q v c 



APPENDIX C 


POSITION OF CURVED REENTRIES 


In figure 6, iq define line U constrained to pass through the center of 

projection of camera A and some point on the curved reentry. These direction cosines 
are obtained by evaluating the polynomial fit through the reentry image at some point 
xj,yi and using equations (10) and (12). The equation of the line U relative to the 
astronomical equatorial system X,Y,Z can be written as 


X_ = Y_ = Zl 
M Mi Vi 


(29) 


Let u be some general point on line U. The coordinates of u are X u , [ ~i)X u , 

/n\ ' ll 

I ^7jX u . Line Q constrained to pass through the center of projection of camera B 

(chopping camera) and some point on the reentry corresponding to some shutter break 
on photographic plate B is defined by direction cosines The equation of 

line Q becomes 


X - X 


AB 


Y - Y AB _ Z - Z AB 


Some general point q on line Q has coordinates Xq, 

V\ 

(Xq-X AB )^-+ The distance between the two points 


( X q- X AB)^ +Y AB- 
u and q is given by 


D = 


( X u- X q) + 





Expressions for X u and 
which give a minimum value of 


Xq, in terms of the direction cosines of the two lines, 
D, are determined from 


and 


9D 

3X. 


= 0 


u 
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Solving the two resulting equations for X u and X yields 

4 


where 


X u = 




X, 


_ X j( B i A ij ~ C j) 


1 - A 


+ X 


ij 


AB 




B i “ x i x AB + ^i Y AB + ^i Z AB 
C j = X j X AB + M j Y AB + v j Z AB 


and 


A ij = *i x j + MiMj + 


( 30 ) 


Substitution of X u and Xq as given by equations (30) into the equation for D gives 
rise to an expression for the minimum distance between the two lines defined by direc- 
tion cosines and Xj,/ij,i/j. The resulting expression for D is a function of 

the direction cosines and the relative location of the two cameras. 


The set of direction cosines which gives D a minimum value is found 

by using the numerical scheme discussed in the text. The X,Y,Z coordinates of the 
jth shutter break, relative to camera A, are then determined from equations (29) 
and (30). Hence, 



X 


i ( B i ~ c j A ij) 

1 ' N 2 



m x 

Xi X Aj 


Z 


AJ" 



Aj 


and the slant range from station A to the intersection point is 


R Aj - ( X Aj 2 + V + Z Aj 2 


1/2 Bi - ChAj 




1 - A - 2 

1 Ay 
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Figure 7.- Axis systems associated with the point-by-point triangulation method. 



z 



Figure 8.- Relative station coordinates at instant of reentry. 
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